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We use trajectory calculations to successfully explain two-photon 
"ghost" diffraction, a phenomenon previously explained via quantum 
mechanical entanglement. The diffraction patterns are accumulated one 
photon pair at a time. The calculations are based on initial correlation 
of the trajectories in the crystal source and a trajectory-wave ordering 
interaction with a variant generator inherent in its structure. Details 
are presented in comparison with ordinary diffraction calculated with 
the same trajectory model. 



I. INTRODUCTION 



In the "ghost" diffraction experiment of [|l|, photons correlated in a down conver- 
sion process are separated and then counted in coincidence after one photon passes 
through a slit mask. Strange as it may seem, the slit diffraction pattern is observed 
on the detection screen of the other photon if one only counts those pairs for which 
the first photon is detected by a small-aperature fixed detector. For this EPR 0] 
type experiment, the data was interpreted in in terms of entangled states. A sec- 
ond and similar experiment exhibited optical imaging via two photon coincidence 
counting. This experiment was likewise analyzed in terms of entangled states. In an 
earlier conference report 0] we have shown with detailed numerical calculations that 
"ghost" diffraction patterns can be explained by trajectory correlations assuming an 
initial position and velocity relationship for the particles. 
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In the present work, we give some detailed calculations of "ghost" diffraction 
patterns and examine trajectory details in comparison with those involved in deter- 
ministic explanations of ordinary diffraction. The trajectory calculations are carried 
out numerically using the trajectory- wave ordering interaction of The latter is 
a set of second order differential equations for trajectories which invoke, via solution, 
the integral method of generating random variants. We begin with a brief description 
of this particle-wave interaction. 



II. TRAJECTORY- WAVE ORDERING INTERACTION (TWOI) 

The TWOI studied earlier by the author, made use of a general line parame- 
ter PJ^. The following equations represent a version with time as the independent 
variable. 

dxi dui , . , , > 

— ; ■;- = —{Xi — Ui)F (1) 

dt dt ^ ' ^ ' 

In these equation the Xi represent the particle coordinates, and Xi = dxi/dt. The Ui 
are components of a field ray-velocity and the function F is a real positive function 
of field components with units of inverse time. The field affects the trajectories 
only through the Ui and function F. The TWOI equations can be implemented with 
different choices for the Ui and function F, facilitating the study of different trajectory 
models. Specific expressions for Ui and F used in the diffraction calculations here are 
discussed in the following section. Modifications to impose metric constraints are 
discussed in P,^]. 

The above set of trajectory equations has in general two important properties. 
First, the Xi approach the limit. Second, in this limit (xj ^ Uj) , ( named 

the " attractor limit" in ^ ) the large number distribution of the particles approaches 
the function F provided the initial particle phase-space values are chosen uniformly 
randomly over a large but not unique volume. The attractor limit is approached 
closer and closer as time increases. For practical calculations, one must set some 
criteria for "sufficiently close". These feature can be more easily seen with a simple 
one-dimensional case with Ui = . We directly integrate the following equation. 



This gives the following integral relations 



x{t) = x{0) - / Fx{t')dt' (3) 







x{t) = x{0) - / Fdx' (4) 



If we numerically solve @ until the attractor limit x(t) — » is sufficiently satisfied, 
we obtain the following relation. 



X 



Fdx' ^ x(0) (5) 

'xo 

Repeating this process many times with {x{0),x{0)) chosen uniformly randomly over 
a sufficient volume of phase space, we generate a random variant specified by F. The 
quantity Fdx is proportional to the number of particles that reach the attractor limit 
x{t) — s> between x and x + dx. It is a technique that indirectly invokes the integral 
method to generate in the attractor limit a random variant specified hj F . It re- 
arranges an initial distribution to match F. This interaction is viewed as one possible 
mechanism by which nature could bring about the observed relation between particle 
distributions and wave intensities. 

To solve the trajectory equations we use both fourth and fifth order Runge-Kutta 
methods [§,§1,01 with adaptive step size techniques [jTU|,|rT|. These two methods are 
used as checks against each other. The adaptive step technique saves two to four 
orders of magnitude on net computation time. The system of second-order differential 
equations is reduced to a system of first-order differential equations following the 
method of Butcher . Numerical details and a test case are given in the appendix 



of [0. Figure |I| presents a histogram of particle positions in the attractor limit 
generated using (||) with x(0) = —7 , F = 0.3 exp(— .3|x|) cos^(2x) , and x{0) chosen 
uniformly randomly from the range (0, 1.2). The "sufficiently close" to the attractor 
limit criteria used was \x{t)\ = 10~^ max |i;(0)|. The solid line in Fig. |l| gives the 
profile of F scaled to match the particle histogram. This distribution can also be 
generated with a wide variety of choices of the initial phase space volume. Use of 
random or a variety of amplitudes in F (instead of 0.3 ) together with a sufficiently 
large initial phase-space volume will still generate the same shape. The noise in the 
distribution is due to initializing with a random number generator, the "sufficiently 
close" criteria and numerical approximations. 



The second property of the TWOI equations is the ahgnment feature of the particle 
velocities with the wave ray velocity. From (0) we can see that the rate that the 
attractor limit is approached depends on the initial velocity difference (xj — ui) and 
upon F. In Fig. Q we give some representative trajectories with particles initialized 
at Z = lOA beyond a slit of width lOA. The particles are initialized in directions 
at 20" intervals. One can see that the particles initially moving backwards can be 
turned around and approximately aligned up with the wave ray velocity in less than 
a wavelength. 

The TWOI approach has some similiar features to, but some conceptional differ- 
ences from, the deterministic approach initiated by Madelung Ijl^ and de Broglie |15 



and re-considered later by Bohm 0,0] and others [|18[ for the Schrodinger equation. 
In that work, the particle velocity x is identified with the wave ray velocity so that 
the particles follow the wave flow tangent lines. We will refer to this case { x = u) 
as the attached case. This attached case is approached in the attractor limit for the 
TWOI. 

Critics of attached models quickly point out that the particle distribution must 
initially match the quadratic wave density, and that there is no mechanism to bring 
about the attached condition [x = u). If initial conditions are met, the continuity 
equation guarantees that the distributions will continue to match. However, if tangent 
lines are initiated with a density that does not match the quadratic wave density, they 
will continue to differ. Interestingly, it was Bohm [|T6| who suggested using a force 
that depended on the difference { x — u ) to bring about the attached 
limit. The TWOI can bring about this case as a limit, as well as generate in the 
particle distribution the wave intensity pattern if the latter is represented by F in the 
TWOI equations. The TWOI is a somewhat more general approach in that it does not 
depend upon a continuity equation, and F can represent any field intensity. Likewise, 
the u can be chosen in different ways, e.g., as a vector parallel to the Poynting vector 
as in [O] , or as in the interesting work of Floyd pO| . 



III. SLIT DIFFRACTION 



Slit diffraction has been observed for many types of particles, including photons, 
neutrons and atoms. Can slit diffraction be explained via calculating particle tra- 



jectories, one particle at a time? One adamant view on this question is reflected in 
the well-known comment of Richard Feynman " a phenomenon which is impossible, 
absolutely impossible, to explain in any classical way, and which has in it the heart 
of quantum mechanics " By treating the wave as a field and using the TWOI, 
the author has been able to obtain slit diffraction patterns, one particle at a time. 
Here, we present a single, a double and a triple slit diffraction pattern in Figures |^, 
H, and H respectively. These diffraction patterns are calculated using the following 
scalar wave equation V^\I' — k^"^ = and standard Rayleigh-Sommerfield diffraction 



with Neumann boundary conditions together with the following equation for u. 

A ^ 

u = -z_(^tv^ - (V*^)*)/(^^^) (6) 

For light diffraction, we use A = c/k where c is the speed of light. With (j^) we 
are taking u proportional to the gradient of the phase of the complex field, as in 
the work of Madelung fl^ and de Broglie The scalar field, together with its 

first and second partial derivatives needed in (|1|), are calculated numerically for each 
point along each trajectory. Because an adjustable step technique is used, the exact 
position where the trajectory crosses a screen is determined by using a fifth order 
Lagrange interpolation polynomial . 

In each figure, the wave is calculated assuming a point source 500A in front of the 
slits and a detector screen at lOOOA beyond the slits. For the single slit case shown in 
Fig. 0, the slit width is 8A. In the double slit case shown in Fig. |^ , the slit width is 5A 
and the center-to-center slit separation is 25A . In the three slit case shown in Fig. ^ , 
the slit width is 5A , and the center-to-center slit separation is lOA. The solid lines in 
these figures represents the wave intensity scaled to match the histograms. From these 
figures, it is clear that slit diffraction patterns can be reproduced via deterministic 
trajectory calculations using the TWOI and treating the wave as a field. 

It is instructive to view some typical trajectories. Figure |^ shows some typical 
trajectories for a three slit case with slit width of 5A and slit center to center sepa- 
ration of 7A. These trajectories are near the attractor limit at the slits, so that the 
trajectories closely represent the wave flow tangent lines. Notice that these trajec- 
tories are not straight lines and that the trajectories do not cross each other. These 
feature are characteristic of tangent line trajectories for waves involving the operator 



IV. TWO-PHOTON GHOST DIFFRACTION 



Two-photon correlations involving slit diffraction have been measured using pho- 
ton pairs made in a parametric down conversion process [jl| . Figure ^ illustrates the 
experimental setup. The "unfolded" sketch is illustrated in Fig. These figures 
and notation follow Fig. 5 of Ref. with some changes that feature the trajectory 
calculations here. The two polarization components are separated by a polarizing 
beam splitter into the "signal" beam and the "idler" beam as shown. The signal 
beam passes through a slit pattern onto a screen about a meter behind the slits. 
No slit pattern is placed in the idler beam. A small aperture detector Di is placed 
at a fixed position on the signal screen and a small aperture detector D2 is used to 
scan the idler screen. Coincident two-photon measurements are made for signal-idler 
photon pairs if a signal photon enters the fixed detector Di. After many pairs have 
been measured, the accumulated distribution on the idler screen clearly exhibits a 
slit diffraction pattern Q]. This is remarkable because no slit mask exist in the idler 
beam. 

The author considers this to be one of the most intriguing experimental results 
of modern physics. It would appear that what happens to the signal photon is some 
how "known" by the idler photon. A two-photon entangled state explanation of 
the data distributions was given in Can this type of experiment be explained 
with a deterministic interactive particle-field model? To investigate this problem, 
we consider a very simple model whose features are best illustrated by the unfolded 
version sketch in Fig. ^ b. In this figure, as in , the distance Z2 is the distance 
from the slits back through the beam splitter BS to the crystal and then back along 
the idler path to the idler screen. The distance Zq is the distance from the crystal to 
the beam splitter and on to the slits. Zi is the distance from the slits to the fixed 
detector Di. 

We describe this system using our interacting particle-field model. The o-ray 
and e-ray polarization components satisfy the wave equations V^-Eo — k'^Eo = and 
V^-Ee — k'^Ee = 0, and for both components we assume spherical wave solutions with 
sources centered at some point S in the crystal. The e-ray is diffracted by the slits 
and the o-ray continues to the idler screen. We likewise assume that together with 
the spherical wave, two particles are initialized at point S in the crystal with equal 



and opposite velocity vectors For the e-ray we calculate the diffracted wave as in 
the slit diffraction described above, but here the source point S is generally not on 
the optical axis. 

To calculate the particle trajectories we use Eqs. (|l]) and with A = c/k, and 
\1/ = -Eo or \1/ = i^e- Our model crystal source is very thin in the Z direction and for 
each pair of particles, we choose the Y position of the source point S in the crystal 
uniformly randomly from a range of {—W,W). The initial angular direction of the 
signal particle is likewise chosen uniformly randomly. However, the initial velocity 
direction of the idler particle is always fixed as the negative of the initial velocity 
direction of the signal particle (as viewed in the unfolded picture of Fig. |^). This 
and the same initial position represent the initial trajectory coupling that gives rise 
to the "ghost" diffraction in this model. We first calculate the trajectory of the signal 
particle, and if it terminates in the aperture of detector Di, we calculate the trajectory 
of the idler particle and record where it hits the idler screen. With this coincidence 
restriction, the accumulated histogram of idler particles on the idler screen is the 
"ghost" pattern in the language of An initial report by the author on this type of 
two-photon "ghost" diffraction calculations for a single and double slit case appears 
in a recent conference proceedings Q]. 

We show in Fig. ^ and Fig. ^ the "ghost" diffraction patterns calculated with this 
simple model for cases of one and three slits respectively. For both cases we used 
Z2 = 7000A, Zq = lOOOA, and Zi = 500A. For the single slit case, shown in Fig. || we 
used a slit width of 15A, a crystal half width oiW = 155A , and a 6A aperture for the 
signal detector Di. For the three slit case shown in Fig. we used a slit width of 5 A, 
a slit center-to-center separation of 7A, a crystal half width oiW = 215A, and a lOA 
aperture for the detector Di. It was necessary to calculate many signal trajectories 
because very few terminated on the narrow detector aperture of Di. 

The distinct recognizable diffraction patterns exhibited in Fig. ^ and Fig. ^ clearly 
indicate that two-photon "ghost" diffraction can be explained via deterministic tra- 
jectory calculations using the interacting particle-field model. This two-photon corre- 
lation phenomenon is a result of pattern selection via correlated trajectory selection. 
The trajectory correlations are entirely due to initial conditions. In Fig. ^ and Fig. ^ 
the solid line represents the diffraction pattern on the idler screen that results as if 



one had a point source at the signal detector on the signal screen. The common 
features of the solid line and the histogram are due to the common geometry in the 
almost Fraunhofer limit. However, the patterns differ slightly from each other. This 
difference can be reasonably explained by the fact that the signal detector Di has 
a finite size. It is conceptually incorrect to think that the "ghost" diffraction is the 
result of a "collapsed probability" signal at Di propagating backwards to the idler 
screen. 

In the conventional slit diffraction of Fig. ^, Fig. ^ and Fig. ^ the wave used 
for each particle had a common point source. By sharp contrast, each particle in 
the "ghost" diffraction generally is influenced by a different wave since the source 
positions in the crystal have a wide range {—W, W) . A random selection of signal 
trajectories just beyond the slits is shown in Fig. |TD| for the three slit case. One can 
see that the trajectories cross each other and appear to have no definite pattern. In 
fact, if one accumulates a histogram on the signal screen, one finds a broad smear 
with no definite pattern. This is expected because each wave generally starts from 
a different source points in the crystal. However, it is in this very situation that we 
obtain the "ghost" diffraction pattern on the idler screen for that selection of signal 
trajectories that terminate in the fixed detector Di. This feature is consistent with the 
observations in IQ, in which the "ghost" patterns were obtained using rather divergent 
beams. By shifting the position of the detector Di we obtain a shifted "ghost" pattern 
on the idler screen. This corresponds to selecting an entirely different set of signal- 
idler trajectory pairs. The authors in [||] mentions that this interesting feature was 
observed in the experiment. 



V. SUMMARY 



We have made "ghost" diffraction calculations using a trajectory model with two 
key features. The first feature is the use of the TWOI equations to describe the 
trajectories. Through the TWOI, the diffracted field influences the distribution of 
the trajectory particles, one trajectory at a time. As a physical model both particle 
and field must be taken together. Here, the TWOI used to calculate the particle- 
field interaction is much closer to that of classical physics in which the particle is 
not "attached" to the field. This key freedom, not found in the "attached" models 



initiated by Madlung and de Broglie many years ago, facilitates the random 
variant generator inherent in the structure of the TWOI to generate the wave intensity 
patterns in the ensemble distribution. 

The second key feature is the initial velocity correlation of the photons in the 
down conversion crystal, as well as the same initial position. As the calculations 
clearly demonstrate, this initial coupling together with the TWOI is sufficient to give 
rise to the "Ghost" diffraction patterns. This is a simple classical trajectory model 
of this photon-photon experiment. Its success clearly lends support to a realistic and 
deterministic description of nature. 
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FIG. 1. A one dimensional histogram of particle count versus the X position in the 
attractor limit. The solid line is the shape of the function F scaled to approximately 
match the histogram. 
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FIG. 2. Typical trajectories illustrating how the TWOI can bring about the 
attr actor limit (xj — > Uj). The trajectories are initialized at the same point in 
directions at 20^^ intervals. 
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FIG. 3. Histogram of particle count versus the Y position on the detector screen 
for single slit diffraction. The trajectories of the particles were calculated one 
at a time. The solid line represents the wave intensity curve scaled to approximately match 
the histogram. 
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FIG. 4. Histogram of particle count versus the Y position on the detector screen 

for double slit diffraction. The trajectories of the particles were calculated one at 

a time. The solid line represents the wave intensity curve scaled to approximately match 

the histogram. 
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FIG. 5. Histogram of particle count versus the Y position on the detector screen 
for three-sht diffraction. The trajectories of the particles were calculated one at 
a time. The solid line represents the wave intensity curve scaled to approximately match 
the histogram. 
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FIG. 6. A selection of typical trajectories in the attractor limit in the very- 
near field region for three-slit diffraction. The dots on the lines represent the 
adjustable step positions. The trajectories are not straight lines and do not 
cross each other. 





FIG. 7. Schematic (a) illustrates the experimental setup and calculation model for 
"ghost" diffraction, and schematic (b) illustrates the unfolded view following Fig. 5 of Q. 
Detectors Dl and D2 measure the signal and idler photons in coincidence. Dl is fixed 
and D2 is scanned. The diffraction pattern is observed on the idler screen whereas the 
signal screen is behind the slits. 
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FIG. 8. A random selection of signal trajectories just beyond the three-slit mask. 
Even though the particle is each case is near the attractor limit, the trajectories can 
cross each other because each is influenced by a wave from a different point source 
location in the crystal. The "ghost" diffraction pattern arises from those idler 
trajectories paired with just the very small fraction of signal trajectories that 
terminate on the small detector Di. 
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FIG. 9. A single slit "ghost" diffraction pattern. Histogram of photon pair count 
versus Y position on the idler screen. The correlated signal and idler trajectory 
pairs arc calculated one at a time. The solid line represents a diffraction 
curve calculated on the idler screen as though one had a point source 
at the signal detector. 
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FIG. 10. A three slit "ghost" diffraction pattern. Histogram of photon pair 
count versus Y position on the idler screen. The correlated signal and idler 
trajectory pairs are calculated one at a time. The solid line represents a diffraction 
curve calculated on the idler screen as though one had a point 

source at the signal detector. 




